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The transverse momenta of hadrons in central nucleus- 
nucleus collisions are evaluated in a boost invariant hydro- 
dynamics with transverse expansion. Quark gluon plasma is 
assumed to be formed in the initial state which expands and 
cools via a first order phase transition to a rich hadronic mat- 
ter and ultimately undergoes a freeze-out. The average trans- 
verse momentum of pions, kaons, and protons is estimated for 
a wide range of multiplicity densities and transverse sizes of 
the system. For a given system it is found to scale with the 
square-root of the particle rapidity density per unit transverse 
area, and consistent with the corresponding values seen in pp 
experiments at 1800 GeV, suggesting a universal behaviour. 
The average transverse momentum shows only an approxi- 
mate scaling with multiplicity density per nucleon which is at 
variance with the pp data. 

PACS numbers: 12.38M 



The last two decades have witnessed an active pursuit 
of relativistic heavy ion collisions in order to create a 
quark-gluon plasma, first 'seen' a few micro-seconds af- 
ter the big-bang. The recently concluded experiments at 
the CERN SPS and the ongoing experiments at the Rel- 
ativistic Heavy Ion Collider at Brookhaven do seem to 
suggest a very definite possibility of the confirmation of 
the quark-hadron phase transition. 

It is known for quite some time jjj that the transverse 
momenta of particles produced in a high multiplicity en- 
vironment in almost all models of particle production, 
roughly scale with the particle rapidity density (dN/dy) 
per unit transverse area (ttR%), 



assumption is made about the applicability of Bjorken 
hydrodynamics H we get, 
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This follows from dimensional arguments as v is the only 
scale in the problem. This has been noted in minijet mod- 
els H , string models Q , and more recently suggested in 
partonic saturation descriptions like colour gluon conden- 
sate model for heavy ion collisions Q] . 

Let us approach this problem from another angle. Let 
us assume that in central nucleus-nucleus collisions at 
relativistic energies a thermally and chemically equili- 
brated quark gluon plasma is formed at an initial time 
tq at an initial temperature Tq. Once v is known and an 
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where a = 42.257r 2 /90 for a plasma of massless quarks 
(u, d, and s) and gluons and we have put the number of 
flavours as « 2.5 to account for the mass of the strange 
quarks. Do the hadronic spectra resulting from the hy- 
drodynamic expansion of the plasma thus produced show 
some scaling behaviour? 

Assuming that one has (E) = 3T for massless parti- 
cles, we may take Tq = 1/3T , from considerations of the 
uncertainty relation. Thus we see that the initial tem- 
perature scales with u 1 / 2 . We may also add that it was 
suggested some-time ago that HQ the S+Au and Pb+Pb 
collisions at the CERN SPS present an interesting exam- 
ple as they lead to systems with different transverse sizes 
but similar initial temperatures due to the near equality 
of the scaling variable v. (We shall see that the differ- 
ences for the hadronic spectra for the two cases are not 
large.) 

In the present work we continue this argument and 
show that the average transverse momenta of the pro- 
duced hadrons, estimated on the basis of a hydrodynam- 
ical expansion with transverse flow scale with v 1 / 2 . The 
scaling is surprisingly similar to the corresponding be- 
haviour for pp data at 1800 GeV, which is sought to be 
understood within a colour gluon condensate model [Q. 
It may be recalled that these data were used to provide 
a suggestion for baryon-flow and quark gluon plasma [|| . 
A straight forward explanation in terms of an increased 
minijet activity which leads to a larger multiplicity has 
also been offered for this ||. 

The similarity of these behaviours is of interest as 
the hydrodynamics calculations start with a partonic 
plasma, which develops into a mixed phase and then into 
a hadronic phase. The spectra for particles having differ- 
ent masses arc affected differently due to the celebrated 
radial flow velocity which develops during the expansion. 
More-over several details like the equation of state, the 
critical temperature, the freeze-out temperature, and a 
model of hadronization etc. affect the overall evolution 
of the system. 

This would seem to suggest that the scaling of the par- 
ticle transverse momenta may be 'seeded' in the initial 
conditions, and that this behaviour is broadly retained 
by the hydrodynamic expansion. We insist that radial 
flow is essential to obtain this result, as in absence of the 
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radial flow, the (j>t) is uniquely determined by the freeze- 
out temperature and is independent of multiplicity Jll]] 
and the size of the system. 

The method of solution of the hydrodynamic equations 
for the expansion and cooling of the quark-gluon plasma 
has been discussed extensively in literature. We only 
briefly mention the salient features. As discussed earlier, 
we assume the plasma to be in a thermally and chemically 
equilibrated state of quark gluon plasma, which expands, 
cools, undergoes a phase transition to a hot hadronic 
matter at T — Tq and when the hadronic density is too 
low to cause further scattering among the hadrons it un- 
dergoes a freeze-out. The plasma is assumed to undergo 
a longitudinally boost-invariant and an azimuthally sym- 
metric transverse expansion. 

The initial energy density profile is assumed to fol- 
low the so-called 'wounded- nucleon' distribution flipf . We 
further assume that the phase transition takes place at 
T = 180 MeV and the freeze-out takes place at 120 
MeV. We use a hadronic equation of state consisting of 
all hadrons and resonances from the particle data table 
which have a mass less then 2.5 GeV Oj. In all the 
calculations reported here the initial time tq is taken as 
1 fm/c. It is well known that as long as TqTq is kept 
fixed, a smaller initial time makes only a marginal differ- 
ence to the flow. The relevant hydrodynamic equations 
are solved using the procedure fl2|l discussed earlier and 
spectra for hadrons estimated ||11[] using the Cooper- Fry 
formulation. 

In Fig.l we show our results for several systems |^| hav- 
ing an identical particle rapidity density per unit trans- 
verse area (v). We see that the average transverse mo- 
menta for pions, kaons, and protons for all the cases are 
similar, though they decrease slightly with the increas- 
ing transverse size. This decrease is about 4% for pions 
and about 7% for protons. Even though marginal, this 
is interesting as the radial flow introduces an additional 
scale in the problem, the size Rt- At a time t the re gion 
beyond Rt — c s t is affected by the transverse flow [ O . 
Thus a smaller system is more strongly affected by the 
flow, even if the initial temperature is similar and thus 
build-up Jul of the radial flow velocity given by 
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will be affected. In the above, P is the pressure and 
e is the energy density so that the quantity under the 
integral sign gives the speed of sound which passes though 
a minimum during the mixed phase. 

A demonstration of this effect is seen in Fig. 2, where 
we have plotted the emission of pions per unit transverse 
area for the S+Au (WA80 experiment |l5ll) and Pb+Pb 
(WA98 experiment |l6|]) systems which have the same 
v but different transverse sizes (Rt)- We add that the 
hydrodynamic treatment used here gives a perfect de- 
scription to the hadronic spectra for the Pb + Pb system 
at SPS energy, which correspond to this value of v p0|. 



We note that the inverse slope of the smaller system is 
marginally larger, as indeed Fig.l implies. The slight 
difference in the radial flow which may develop in sys- 
tems of varying sizes can be magnified by looking at the 
transverse momenta of heavier particles. This is seen in 
Fig. 3 where we now plot the emission of protons per unit 
transverse area for the S+Au (NA35 @) and Pb+Pb 
(NA44 fig] ) systems having the same v. A much larger 
freeze-out temperature ( with very different values for pi- 
ons and protons) would be needed to understand these 
spectra if the flow is not included. 

In Fig. 4 we show the average transverse momenta of pi- 
ons, kaons, and protons for central collision of lead nuclei 
for a range of multiplicity densities spanning the region 
likely to be covered at the RHIC and LHC. It is seen that 
the average transverse momentum rises with the square- 
root of the scaling variable v. so that, the pion, kaon, 
and proton spectra determined from the hydrodynamics 
are well described by: 



(p T ) = 0.326 + 0.052^ ; tt , 
(p T ) = 0.423 + 0.087V^ ; K , 
(p T ) = 0.530 + 0.137^ ; P ■ 



(5) 



We also find (see Fig. 5) that the the hadronic-spectra 
scale as 



dN h a 
d 2 p T dy (p T ) 



Pt 
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The preliminary data for the transverse momenta are 
now available from the PHENIX experiment |l{| at 
^snn 130 GeV. We see from the Fig. 6 that these data 
conform to the hydrodynamic behaviour seen in this 
work. 

In Fig. 7 we compare the hydrodynamics results for 
Pb + Pb system with the experimental results for the 
pp collision at 1800 GeV |pf|. Their agreement points 
to a universal scaling behaviour fot the transverse mo- 
mentum distributions. Considering that several models 
of particle production relate the increase in (pt) > with 
the multiplicity to the initial conditions, we have plotted 
the (pt) of the hadrons against the initial temperature 
obtained and used in the hydrodynamic calculations here 
and we see (Fig. 8) that the increase is directly related to 
initial temperature. 

It has been suggested |2l[] that hydrodynamic flow at 
the end of the three dimensional evolution of the plasma 
is uniquely determined by a dimensionlcss parameter ui: 



s t /s c Rt 



(7) 



where sq is initial entropy density, tq is initial time, s c 
is the entropy density at the beginning of the phase- 
transition, and Rt is the transverse size. Recalling that 
s ~ T 3 , one can easily prove that ui — t c /Rt, where 
t c is the proper time when the phase transition begins, 
provided the QGP cools according to a boost-invariant 
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longitudinal expansion. Obviously if ui is small, trans- 
verse flow effects should be small, as the QGP phase is 
over well before Rt / c s , when the rare- faction wave cov- 
ers the distance Rt and agitates the entire fluid. Using 
Eq.(||), we can also see that lu ~ (dN/dy)/A. We now 
investigate this scaling in terms of a variable 



/i ~ A dy 



(8) 



in the spirit of a similar study by Ruuskanen ]2^] . 

In Fig. 9 we have plotted the average transverse mo- 
menta of hadrons obtained for S+S, Zr+Zr and Pb+Pb 
systems for a large range of values of the scaling variable 
li. We see that the professed scaling is only approxi- 
mately satisfied, and becomes worse with increasing mass 
of the hadron at large [i. This could be due to the fact 
that for such values the development of the flow after the 
conversion to hadrons starts, can not be ignored. We also 
note that the this scaling is not universal as the the cor- 
responding results for the pp data show a very different 
behaviour. 

Before concluding, let us discuss the uncertainties in 
our estimates. The hydrodynamics calculations have sev- 
eral inputs which can slightly alter the extent of the devi- 
ations of the scaling behaviour observered here. We have 
used a freeze-out temperature of 120 MeV. Decreasing Tp 
increases the (pt) slightly, the largest increase being for 
protons. For decreasing Tp to 100 MeV increases 

the (pt) for pions from 0.44 GeV/c to 0.46 GeV/c, while 
for protons the corresponding increase is from 0.83 GeV/c 
to 1.02 GeV/c. For v = 50, the (pt) for pions increases 
from 0.70 GeV/c to 0.72 GeV/c when the freeze-out tem- 
perature is decreased to 100 MeV, while the increase for 
protons is from 1.5 GeV/c to 1.7 GeV/c. The critical 
temperature used affects the {pt) too, in a very interest- 
ing manner. Increasing the Tc increases the degrees of 
freedom in the hadronic matter (as higher resonances get 
excited) and the phase-transition is completed quickly. 
Thus the duration of the low-pressure gradient associ- 
ated with the mixed phase is reduced and (pt) increases 
slightly, and the opposite happens when a lower value of 
Tc is used. We have verified that the change in (pt) for 
protons is less than 10% when Tc is changed by ± 20 
MeV. For pions the results are much less altered. It is 
easy to argue that if a much simpler equation of state for 
hadrons is used (say having only pions) then the mixed 
phase will live much longer and (pt) will decrease. 

It is also useful to remember that the pt broadening 
and the so-called jet-quenching should also affect results 
discussed in Fig. 3, e.g. It is hard to decide the relative 
worth of hydrodynamic and such treatments till precise 
data at higher energies and upto large pt are available. 

The discussion so far has involved central collisions, 
as it is not easy to extend the hydrodynamic treatment 
to non-central (non-zero impact parameter) collisions. A 
rough estimate could still be made if we ignore the az- 
imuthal asymmetry of the region of overlap, approximat- 
ing it to a circle having a radius R ~ 1.2(Vp ar ^/2) 1 / 3 , 



where iVp ar t is the number of participants. This would 
then help assign a value of v even to non-central colli- 
sions, as has been done in Ref. [Q. Recall that the NA49 
data for Pb + Pb collisions treated in this manner were 
found to similar to the E735 data for pp in the above 
work bolstering the arguments for a universal behaviour. 
We eagerly await the data from RHIC experiments to 
verify this aspect. 

The similarity of the behaviour observed here for high 
multiplicity pp and AA results raises one question, which 
is applicable to all the studies involving hydrodynamics; 
viz., are the conditions appropriate for the applicability 
of hydrodynamics realized in these situations? The simi- 
larity of the behaviour for the high multiplicity events in 
pp with those of A A does suggest a plausible argument. 
It indicates that the the initial states of these systems 
are perhaps consisting of partons and the multiple colli- 
sions indeed drive the system to an equilibrium quickly. 
It is known from several studies that high multiplicity 
events in hadron-hadron collisions are also rich in mul- 
tiple partonic collision s |9j | and provide a basis for the 
parton cascade model [E3L 

In brief, we see that hadronic spectra generated from 
boost-invariant hydrodynamics with transverse expan- 
sion show a scaling behaviour in the particle rapidity 
density per unit transverse area. For smaller values of 
the scaling variable i>, the results are quite close to val- 
ues inferred from ISR data, which are sought to be un- 
derstood in terms of a universal behaviour given by the 
colour gluon condensate model. 
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FIG. 1. A model calculation of the average transverse mo- 
mentum of hadrons for several systems which have identical 
particle rapidity density per unit transverse area. 
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FIG. 2. Emission of pions per unit transverse area for 
S+Au (WA80 fn|) and Pb+Pb (WA98 @) systems, which 
have the same multiplicity per unit transverse area but dif- 
ferent transverse sizes. 
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FIG. 3. Ehe emission of of protons per unit transverse area 
for S + Au and Pb + Pb having same v. 
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FIG. 4. A model calculation of the average transverse 
momentum of hadrons for central collision of Pb nuclei as a 
function of the particle rapidity density per unit transverse 
area. The symbols represent the hydrodynamics results while 
the solid curves give the fits described in the text (Eq.H). 



FIG. 5. Scaling of hadronic spectra obtained from hydro- 
dynamics (Eq.^). The two curves differ by a factor which 
is equal to the ratio of the square of the average transverse 
momenta (see text). 
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FIG. 6. Hadronic spectra from PHENIX experiment |l9| . 
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FIG. 9. The average transverse momenta for pions, kaons, 
and protons for central collisions of S + S, Zr + Zr, and 
Pb + Pb systems as a function of multiplicity per nucleon. 
The corresponding values for the pp data Eol are also given. 
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FIG. 8. The average transverse momenta attained as a 
function of the initial temperature. 



FIG. 7. A comparison of scaling suggested by hydrody- 
namics (solid curves) and the corresponding data for pp scat- 
tering at 1800 GeV Q. 
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